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ABSTRACT 

We present optical, X-ray and gamma-ray observations of GRB 111209A, observed at a redshift 
of z = 0.677. We show that this event was active in its prompt phase for about 25000 seconds, 
making it the longest burst ever observed. This rare event could have been detected up to z ~ 1.4 in 
gamma-rays. Compared to other long GRBs, GRB 111209A is a clear outlier in the energy-fluence and 
duration plane. The high-energy prompt emission shows no sign of a strong black body component, 
the signature of a tidal disruption event, or a supernova shock breakout. Given the extreme longevity 
of this event, and lack of any significant observed supernova signature, we propose that GRB 111209A 
resulted from the core-collapse of a low metallicity blue super giant star. This scenario is favoured 
because of the necessity to supply enough mass to the central engine over a duration of thousands 
of seconds. Hence, we suggest that GRB 111209A could have more in common with population III 
stellar explosions, rather than those associated with normal long gamma ray bursts. 
Subject headings: Gamma-ray burst: individual (GRB111209A) 
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1. INTRODUCTION 

Gamma-ray bursts (GRBs) are detected as brief flashes 
of hig h-energy photons typica lly lasting some tens of sec- 
onds (Klcbcsa del et all Il973l ). First discovered in the 
1960s, these cosmological sources are the most ener- 
getic explosion s in the Universe since the Big Bang (see 
lMeszaroa l2006. for a review). There are two main classes 
of GRBs with well se parated properties and behaviors: 
long and short GRBs (|Kouveliotou et aPl 19931 ). Because 
of their extreme distances and diverse emission character- 
istics, it has been difficult to isolate a single progenitor re- 
sponsible for all long duration GRBs. However, the asso- 
ciati on of rare type Ic sup ernovae with several long GRBs 
(e.g. iHjorth et al. l [2003) suggests that these GRBs indi- 
cate the formation of compact objects following stellar 
collapse. This hypothesis, based on the Collapsar model 
(|Wooslevlll993l) . suggests that the proge nitor of long du- 
ratio n GRBs are Wolf-Rayet stars (e.g. FChcvalie r fc Lil 
1999). In addition, other astrophysical phenomena can 
mimic GRBs, such as th e magnetic reconnection onto a 
neutron star (|Usovlll992l ) or the tidal disruption of a mi- 
nor body by a compact object (jBoer et al.lll989h . These 
events must occur at closer proximities than classical long 
GRBs as their available energy is lower. 

We here present an analysis of the longest burst ever 
observed, GRB 111209A , which had a duration of 
25000 se conds. This burst occ urred at a redshift of z 
= 0.677 (|Vreeswijk e t~aT1l2011|) Our analys is uses the 
observ ations of XMM-New ton (|Jansen et~aL|[200l and 
Swift dGehrels et alj 12001 in the X-ray, Konus-Wind 
(lAptekar et all 119951) in the gamma-ray band, TAROT 
( Klotz et al.ll2009h ~iii the optical, as well as archived and 
public data. In this paper we focus our attention on the 
prompt emission of GRB 111209A to determine the most 
likely astrophysical source of such an ultra- long event. 
An additional paper will discuss the afterglow properties 
in a multi wavelength context (Stratta et al., in prepara- 
tion). 

In the following, we define the classes of bursts with du- 
rations greater than 2 s, 10 3 s, and 10 4 s as long, super- 
long and ultra-long GRBs respectively. We also define 
T90 the time required to radiate 90% of the total energy 
observed in a given band. We assume flat A CDM Uni- 



verse with Ho = 71 km s Mpc 



a. 



0.27, and 



quote all errors at the 90 % confidence level. The outline 
of the paper is as follows. 

In Sections [2] and [3] we present the observations and 
the data reduction respectively. The data analysis is 
described in Section [4] We discuss the duration and 
scarcity of events like GRB 111209A in Section [5] and 
in Section |6] consider plausible scenarios for the progeni- 
tor of GRB 111209A . We present our conclusions in the 
final Section 7. 

2. GRB 111209A 

GRB 111209A was discovered by the Swift satellit e 
at T = 2011:12:09-07:12:08 UT (jHoversten et al.ll201lD . 
This event initiated two triggers of the Burst Alert Tele- 
scope (BAT; triggers number 509336 and 509337). It 
showed an exceptionally long duration, was monitored 
up to T + 1400 s. However, the burst started about 5400 
seconds before To, as shown in the ground data analysis 
of the Konus-Wind instrument, and a re-analysis of the 



BAT data showed GRB 111209A was detectable 150 sec- 
onds before the trigger. Swift did not trigger at the start 
of the event as the burst was not in the field of view of 
the BAT instrument. 

GRB 111209A was also detected by Konus-Wind. In 
order to compare the light curve of this event with other 
bursts, we used the publicly available soft band (21-83 
keV; similar to the Swift-BAT 15-150 keV band) light 
curve. The earliest portion of the gamma-ray signal de- 
tected by Konus-Wind featured a weak broad pulse. The 
gamma-ray signal was then observed as a multi-peaked 
emiss ion up to about To + 10000 s (jGolenetskii et al.l 
[20H . Using the Konus-Wind results, GRB 111209A 
had a fluence of (4.86 ± 0.61) x 10~ 4 erg cm" 



isotropic energy Ej. 



(5.82 ± 0.73) x 10 5J erg 



, an 
and 

an intrinsic peak of the spectrum in the y f u space of 
E p = 520 ± 89 keV (|Golenetskii et al.ll201lD . These val- 
ues are i n agreement to wit hin two sigma with the Amati 
relation ([Amati et al.ll2002l ) which empirically links these 
last two quantities for long GRBs. 

Swift /XRT observations sta rted 425 s after the BAT 
trigger (jHoversten et al.l 120111 ) revealing a bright after- 
glow, observed also by Swift/UVOT in the optical-UV 
bands at RA(J2000)=00 /l 57 m 22.63 s and Dec(J2000)=- 
46d 48' 03.8", with an estimated uncertainty of 0.5 arc- 
sec. The afterglow was also clearly detected by ground 
based instrument s; for example the TAROT-La Silla 
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2011 ). and the GROND robotic telescopes 



20111 ). In addition, we activated a Target 



of Opportunity observation with XMM-Newton, between 
T + 56, 664 s and T + 108, 160 s (see the light curve in 
Fig. [!}. This period covered the end of the prompt phase 
seen in X-ray, a subsequent pla teau phase, and the start 
of the normal afterglow decay (jGendre et al.l[20TTI ). 

3. DATA REDUCTION 
3.1. Swift XRT 

Swift/XRT data were taken as Level 1 event files from 
the Swift Data Center archive at NASA/GSFC and were 
processed with the task xrtpipeline version 0.12.6, ap- 
plying the most updated calibrations. Data acquisition 
started in Windowed Timing (WT) mode and continued 
in this mode up to the fourth orbit (T + 22700 s). XRT 
then switched to the Photon Counting (PC) mode for 
the remaining of the observation (up to 26 days after the 
trigger). Pile- up can be observed in WT mode during 
the first orbit (i.e. up to T + 2000 s), and in PC be- 
tween T + 22700 s and T + 30900 s. As a consequence, 
spectr a and light curves wer e ext racted using th e meth - 
ods of iRomano et all (|2006l ) and iVaughan et alj ([2006), 
excluding a central circle of radius 1 pixels from the PSF 
in WT mode; for the PC mode, the radius was set dy- 
namically betwen 7 and 1 pixels as a function of the flux 
of the afterglow, and the results checked with the PSF 
profile of the XRT. 

All spectra and light curves were then extracted us- 
ing standard filtering and screening criteria (including 
bad columns and bad pixels) , and corrected for telescope 
vignetting and aperture filtering using the ancillary re- 
sponse files generated by the task xrtmkarf. We then 
rcbinncd all files to either match the spectral resolution 
when the data quality was high enough, or to obtain at 
least 30 counts per bin. 
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Figure 1. The light curve of GRB 111209A, presented with a temporal axis which is linear in the left panel (the prompt emission) and 
logarithmic in the right panel (the afterglow emission). X-ray data are in grey (XRT) and black (XMM-Newton) . The Konus-Wind data 
(blue solid line) has been scaled to the X-ray data for comparison. The R-band light curve is constructed from TAROT and Swift-UVOT 
data (color indices are indicated on the plot). 



3.2. XMM-Newton 

We retrieved the Observation Data Files from the 
XMM archive and reprocessed them using the XMM-SAS 
version 12.0.1 and the latest available calibration files as- 
sociated with this version of the SAS. The raw events 
files were processed using the tasks emchain, epchain, 
rgsproc, omichain, omfchain. The PN instrument ob- 
served GRB 111209A in Full Window mode, while the 
MOS cameras setting was on Small Window. Because of 
the brightness of the afterglow, we checked for pile-up in 
the data. We used the task epatplot for this purpose, and 
confirmed that the data were free from pile-up, even if a 
significant fraction of out off time events were observed. 
Because of this latter effect, the background regions have 
been chosen, when possible, far away from the source 
and corrected for the difference of off-axis angle using 
the ARFs. In the special case of the PN instrument, the 
instrumental background is highly variable depending on 
the position on the instrument - this was taken into ac- 
count by using a region where the instrumental back- 
ground was compatible with the position of the source. 

We then checked for flares of solar protons, that can 
induce a huge increase of the background count rate, and 
found some marginal events where the background rate 
was increased by a factor of 5. Due to the brightness 
of the afterglow, these events are not very significant for 
the following analysis: when maximal, the background 
represents around 6% of the total signal at the position 
of the afterglow. 

The event files were then filtered for good photons, us- 
ing FLAG == and Pattern filtering (less than 4 for the 
PN, less than 12 for the MOS). Spectra were extracted 



using the task especget, light curves with the task evse- 
lect. 

4. DATA ANALYSIS 

4.1. Temporal analysis 

GRB 111209A X-ray emission indicates an initial shal- 
low decay phase. Assuming a power law model between 
T + 0.425 ks and T + 10.4 ks, the best fit decay in- 
dex is ap = 0.544 ± 0.003. This phase ends very late in 
comparison with typical X-ray emissions. At the end of 
this phase, a gradual steepening drives the light curve 
to the so called "steep decay" phase, thought to indicate 
the end of the prompt emission and is characterized by 
very fast decrease of the count rate. During this phase 
and the following one, the light curve can be modeled 
with a double broken power law with a steep decay in- 
dex of a s = 4.9 ± 0.2, a plateau with decay index of 
ctf = 0.5 ±0.2, and a final decay index of a a = 1.51±0.08. 

4.2. Spectral analysis 

In the following, all models fitted to the data were 
multiplied by two components in order to take into ac- 
count our Galaxy and the host galaxy photoelectric ab- 
sorption from metals. The equivalent hydrogen column 
density of our Galaxy was set to Njj = 1.48 x 10 20 cm~ 2 
(jKalberla et alJl2005f ). The second absorption compo- 
nent was left free to vary at the redshift of the burst, as- 
suming solar metalicity. We found no temporal variation 
of the intrinsic Nh, with a mean of (2.5±0.4) x 10 21 cm~ 2 . 

We began with an analysis of the prompt spectral emis- 
sion between T and T + 2500 s, seen by Swift /XRT. 
GRB 111209A shows a strong flare at T + 1200 s. Since 
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Figure 2. Swift/XRT WT spectrum. The spectrum is fitted with 
a simple power law (upper panel) or with the addition of a black 
body and a broken power law (lower panel). An improvement fit 
is clearly visible in the bottom panel. 

flares are known to have highly variable spectra, this 
was excluded in order to investigate the underlying spec- 
tral continuum. The spectrum was of exceptionally high 
quality, and was re-binned in order to obtain at least 400 
counts per bin. The background subtracted energy spec- 
trum between To + 425 s and To + f 000 s is inconsistent 
with a simple power law or cut-off power law model. The 
best fitting model was that based on a broken power-law 
added to a blackbody component. The high energy seg- 
ment of the broken power-law is consistent with the BAT 
spect rum at the same time (r = 1.48±0.03: IPafmer et al.l 
20ff), but harder than the time averaged Konus-Wind 
best fit model. We note however that the following spec- 
tra (extracted between 5.0 and 47.8 ks) are softer, consis- 
tent with the time averaged Konus-Wind spectrum. We 
therefore assume that the spectral position of the break 
(i.e. Eb rea k) varies with time, and that this break is due 
to the crossing of a specific frequency within the observed 
band. The black body component is very soft (see Table 
[J), and accounts for about 0.01 % of the total flux in the 
0.5-10.0 keV band. Its luminosity is ~ 1.6 x 10 50 erg s _1 . 
Moreover, this component is not detected at later times. 

The XMM-Ncwton monitoring period covered the 
transition from the steep decay phase to the "normal" 
decay phase, for an interval of 53 000 seconds. Because 
of the greater sensitivity of XMM-Newton in comparison 
with the Swift XRT, we were able to extract the spec- 
tra in six different time intervals (see Tab. [T]). A simple 
power law model was first used, giving an acceptable 
value of the x 2 statistic value. However, the residuals of 



the fit clearly showed systematic excesses at the hard end 
of the spectrum, and thus a double component model was 
preferred. The inclusion of a second component produced 
an F-test probability (the probability that the decrease 
of the x 2 statistic value is due to a spurious effect) of 
the order of 10 -9 . The spectral fit clearly indicates that 
this second component is not dominating, and provides 
a small contribution to the total observed flux. 

We extracted a Spectral Energy Distribution (SED) 
during the XMM-Newto n obse rvation using published 
results from IKann et al.l (|2011D in the optical, in or- 
der to confi rm no thermal emi ssion was present during 
this phase. IKann et al.l (|2011l) already noted that the 
GROND data excluded a thermal component, and indeed 
we have found an optical spectral index of /3 = 1.1 ± 0.2, 
using the optical data alone. Using the X-ray to optical 
data, we found that an absorbed broken power law fitted 
the data well. We note that we can exclude the presence 
of a thermal emission in the far UV or the soft X-ray 
parts of the broad band spectrum. 

As stated previously, this paper focuses on the proper- 
ties of the progenitor of GRB 111209A , and the na- 
ture of this second component will be investigated in 
more detail in a forthcoming paper (Stratta et al., in 
preparation). We note however that this is the first 
time, to date, that such a double non-thermal compo- 
nent has been detected during the steep decay phase. 
Several models can account for it (e.g. see the model of 
IBarniol Duran fc Kumail 12009). The optical light curve 
features a rebrightening during the XMM-Newton obser- 
vation: this component may be linked to the prompt- 
to-afterglow transition, to the optical bump, or even be 
totally unrelated to these components. 

5. DISCUSSION 

5.1. Duration of the prompt emission 

The estimation of the burst duration using classical 
methods (i.e. T90) is complicated by two factors: it is 
well known that the duration depends on the observing 
band (the softer the band, the longer the detected signal); 
and for this event the prompt phase is split into several 
orbits, with gaps in the signal. Significantly, the event 
was recorded by Konus-Wind before To, when the burst 
was not in the field of view of Swift/BAT. Therefore, all 
attempts to use classical methods are flawed. To give an 
estimate of the total duration for comparison with other 
bursts, we thus estimated the burst start and stop epochs 
T start and T end as follows: 

• For T start , we use the start of the detection by 
Konus-Wind, about 5400 seconds before To; 

• For T en d, we use the time of the start of the sharp 
decay, about 20,000 seconds after T . 

The choice of the T en d value is motivated by the fact 
that the sharp decline is usually interpreted as the high 
latitude emission of the prompt phase, and its start 
time corresponds to the tru e end of the prompt phase 
(jKumar fc Panaitescul 12000). Before the steep decay 
phase, the initial XRT and BAT light curves are usu- 
ally well correlated, indicating a common origin. This is 
also the case for GRB 111209A (see Fig. [[]). With these 
assumptions, the total duration of GRB 111209A is at 
least 25,000 seconds. 
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Table 1 

X-ray spectral analysis of GRB 111209A. 



Instrument 
observing 
mode 


Start 
time 

00 


End 
time 
(8) 


Spectral 
Model 


r 


Eq 
(keV) 


Second Spectral 
component 
(eV or energy index) 


V 2 
(dof) 


Swift/WT 

Swift /WT 
Swift/PC 
Swift/PC 
XMM-Newton 


425 

5025 
22700 
28780 
55073 


1000 

7825 
25130 
30880 
61416 


BKN+BB 

PL 
PL 
PL 
2PL 


1.17 ±0.03 
1.52 ±0.06 
1.58 ±0.03 

1.8 ±0.1 

1.9 ±0.1 
0.9 ±0.7 


3.4 ±0.4 


30 ±5 
2.7 ±0.2 


0.98 (181) 

1.13 (590) 
0.87 (50) 
0.78 (32) 
0.89 (151) 


XMM-Newton 


61416 


66416 


2PL 


1 2+ - 7 

L - z -0.9 






0.84 (111) 


XMM-Newton 


66416 


71416 


2PL 


I Q+04 

■"••"-o.s 




o o+0.9 

°- a -o.s 


0.86 (95) 


XMM-Newton 


71416 


76416 


2PL 


i 7+0.5 




o 1+0.8 
d,i -0.4 


1.05 (83) 


XMM-Newton 


76416 


91416 


2PL 


i 7+0.5 

— B-jC 




o q+0.5 
<3 - u — 0.3 


0.89 (228) 


XMM-Newton 


91416 


108599 


2PL 


9 O+0.3 

z -°-0.8 




o 4 +i.3 


1.04 (234) 



We extensively surveyed the Swift, BATSE and Konus- 
wind catalogs to find other ultra-long bursts. The 
most complete and homogeneous catalog o f prompt burst 
prope rties is the BATSE 4B catalog (|Paciesas et al.l 
1999). From this catalog, only few bursts have a du- 
ration larger than 1000 seconds (the super- long ones), 
and none have a duration longer than 10 000 seconds. 
For the case of S wift, we also search ed in the XRT light 
curve repository (jEvans et al.l 120071 ) to obtain the start 
time of the sharp decay for all bursts, and found no burst 
with such a large duration. Thus, GRB 111209A is to 
date the longest confirmed GRB ever observed, and its 
duration is exceptional. 

We note however, for completeness, that one event, 
longer than GRB 111209A, has a debated origin. 
Swift J164449.3+573451 (hereafter J1644+57) has been 
proposed by s everal groups to be a tidal disrup- 
tion event (e.g. iBurrows et alj[201ll: iBloom et al.l 120111: 



iLevan et alj|20lH iTchekhovskov et al.ll2013l) . Neverthe- 
less IQuataert fc Kasenl ([2012D . even if agreeing that the 
most plausible origin was a tidal disruption event, have 
shown theoretically that the explosion of a massive star 
could produce a similar event. 

5.2. Selection effects on ultra-long GRBs 

It is possible that selection effects make only the bright- 
est part of ultra-long burst light curves detectable, thus 
reducing their intrinsic durations and making ultra-long 
GRBs indistinguishable from normal bursts. We investi- 
gated this selection effect in two ways. 

Firstly, we constrained the maximum distance from 
which this burst could be detected. We extracted data 
from the Konus-Wind light curve, re-scaling it to a given 
distance and then added the background of the instru- 
ment (adding statistical fluctuations in the final signal 
to take into account Poisson statistics). We used two 
signal levels: the original one, and a second level an or- 
der of magnitude fainter. The light curve was examined 
for a burst with a rate trigger threshold of five sigma. 
We used two different bin sizes for the light curve (1 
and 2 minutes) , the latter simulating more complex trig- 
ger methods where an image trigger can be considered. 
The simulation was performed 100 times for all redshifts 
between 0.3 and 2.2 (with a step of 0.1). Results are 
presented in Fig. [3l and are quite clear: GRB 111209A 




Redshift 



Figure 3. Detection significance expressed in signal to noise ratio 
units for a GRB 111209A like event as a function of its redshift. 
We present the detection regime of such bursts using light curve 
bin sizes of 1 minute (purple line) and 2 minutes (blue line); and 
a burst 10 times fainter with 2 minutes bin size (black line). The 
dotted red line indicates our detection threshold (5 sigma level), 
showing that even in the optimal case the maximum redshift for 
detection is z = 1.4-1.5. 

cannot be detected at high-z. In fact, this event can be 
detected only within a limiting redshift of z = 1.4. 

Secondly, for positive detections (at z < 1.4), we com- 
puted the expected duration of the burst, using the stan- 
dard GRB procedure to obtain T90. Again, this is done 
for two distinct temporal resolutions and flux thresholds. 
Results are presented in Fig. [4] Surprisingly, we find that 
the duration is larger than 10000 seconds in all cases. 
Thus, an ultra-long burst that could trigger an instru- 
ment would always be detected as an ultra-long GRB. 
The lack of such other bursts in ~ 30 years of archives 
indicates a very rare event in the local Universe. Such 
events may have been more common in the distant Uni- 
verse, but their faintness makes their detection unfeasi- 
ble. 

5.3. The rate of ultra-long GRB in the local Universe 

We estimated the rate of ultra-long GRBs in the lo- 
cal Universe using the properties of GRB 111209A . 
Taking into account the effect of the redshift on 
Swift's sensitivity and k-correction (jCoward etalll20rl 
iHowell fc Cowa rd 2011), this rate density can be inferred 
from the measured peak flux. We used the brightest 
epoch of the GRB 111209A prompt phase, as seen by 
Konus-Wind, scaled to the Swift energy band. This is 
then scaled to the limiting detection flux of Swift to ob- 
tain a maximal redshift for detection for Swift, zn m (see 
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Figure 4. Measured duration of a GRB 111209A like event as a 
function of its redshift. Colors are the same as those in Fig. [3] 
The horizontal red dotted line represent 10 000 seconds. The black 
line is truncated at the point where the errors are larger than the 
measured duration. 



ICoward et al.ll20r2l for details). The corresponding vol- 
ume of detectability is then : 



dV 

dz 



dz 



(1) 



Using the cosmological model provided in the Intro- 
duction, the maximum detectable volume correspond- 
ing to ziim = 1.4 i s V^v = 292 Gpc 3 Follo wing 
IGuetta fe Delia Valid (l200l : ICoward et all (12012ft we 
calculate a rate density: 



R 



1 111 

V'max T Q Tj z 



(2) 



with T the maximum observation-time for the sample (7 
years), ft = 0.17 the sky coverage of the instrument and 
r\ z =0.3 the fraction of Swift bursts with a measured red- 
shift. We obtain 9^g 7 x 10 -3 Gpc -3 yr _1 where the errors 
are the 90% Poisson confidence limits (|Gehrerslll986[ ). 

As only four out of 644 Swift bursts detected between 
the launch of the satellite and March 2012 have prompt 
emission durations greater than some thousands of sec- 
onds, we obtain the final rate for ultra long GRBs to be 
of the order of 4 x ii/644 ~6x 10~ 5 Gpc" 3 yr" 1 . This 
is significantly lo wer than that of the no rmal long GRB 
rate density (see lHowell fe Cowardll20liT ). 

5.4. Comparison with super-long bursts 

We retrieved from the literature the prompt emission 
information available for the 9 super-long bursts (with 
duration larger than 1000 seconds) listed in Table[2J The 
fact that GRB 111209A is harder and more fluent than 
other bursts is another indicator of its uniqueness. In 
Fig. [5] we also present these bursts in the duration- 
fluence plane. GRB 11 1209 A appears as an outlier in 
the duration distribution, but not in the fluence distri- 
bution, where other super-long GRBs can have a similar 
(even if slightly smaller) fluence. We note that supernova 
shock breakouts have a much smaller fluence, weakening 
a possible link between this kind of progenitor and GRB 
111209A . Once corrected for distance (when possible 
- several bursts lacked a good follow-up campaign and 
thus a redshift measurement), this burst also appears to 
be more energetic than the others. 

We also investigated the X-ray afterglow proper- 
ties of these bursts. As shown in previous works 




10' 10* 

T90 (sec) 

Figure 5. Position of GRB 111209A (red triangle) in the GRB 
Flucncc-Duration plane. The green stars are long Swift bursts. 
Blue circles and squares are the recorded super long GRBs, and 
the black dot is J1644+57. Blue squares represent high-energy 
transients attributed to SN shock breakout. Note that all measure- 
ments are converted to the Swift band (15-150 keV) for comparison 
with the Swift catalog - therefore some values will differ from those 
given in Table [2] 
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Figure 6. The absolute flux of GRB 111209A. Following the 
method of Gcnd re et all l|2008al ). we rescale the X-ray afterglow 
light curve of GRB 111209A to a common distance of z = 1. 
T his light curve is repr esented by red squares, while the sample 
of lGendre et al.l ]2008al ) are shown as black dots. 

using infrared, optical and X-ray data, the lu- 
minosity of_aftCTgfows tends to cluster in s evera l 
groups (iBoer fc Gendrel 120001: iGendre fc Bolrl [200 



Liang & Z hang) 120061: iNardini et ail 120061: IGendre et al 
2008allbf ). Ifone includes the X-ray ligh t curve of GRB 



111209A in the sample of IGendre etaLl (|2008al see fig. 
[6|), the event is shown to be marginally under luminous. 
However, this could be interpreted as a normal conse- 
quence of the duration of the prompt phase. In this 
phase, internal shocks dissipate the energy of the fire- 
ball. The longer the prompt phase lasts, the less energy 
remains to be dissipated within the following external 
shock. 

6. NATURE OF THIS EVENT 
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Table 2 

Summary of the prompt properties of all bursts with duration larger than 1000 seconds. 



Quantity 


E p 


Ei so Fluencc Duration 


Distance 


Reference 




(keV) 


(10 51 erg) (erg cm -2 ) (T90, s) 


(redshift) 





GRB 971208 


144 


n.a. 


2.6 x 10" 


-4 


2500 


n.a. 


GRB 020410 


n.a. 


n.a. 


2.8 x 10" 


5 


1550 


n.a. 


GRB 060218 


4.9 


0.062 


1.7 x 10" 


-5 


2100 


0.0331 


GRB 060814B 


341 


n.a. 


2.4 x 10" 


-4 


2944 


n.a. 


GRB 080407 


287 


n.a. 


4.5 x 1Q- 


-4 


2100 


n.a. 


GRB 090417B 


> 150 


> 7.3 


8.2 x 10" 


-6 


> 2130 


0.345 


GRB 091024 


280 


350 


1.1 x 10" 


4 


1200 


1.09 


GRB 100316D 


10-42 


0.049 


5.1 x 10" 


-6 


1300 


0.0591 


GRB 101225A 


38 


n.a. 


2.6 x 10" 


-6 


> 2000 


n.a. 


Swift J1644+57 


n.a. 


n.a. 


8.6 x 10- 


-4 


2160000 


0.354 


GRB 111209A 


520 


580 


4.9 x 10- 


-4 


> 25000 


0.677 



Pal'Shin et al. (2008) 
(Nicastro_et_al 
Campana et al. 
Pal'Shin ct al. 
ParShinetaL 
jicjland^t^l^ 
Golcnctskii et al 
J!tarlingj5t_al 




Campana et al. (2011) ; Thone et al. (2011) 
(Burrows et al. 2011) 
This work 



In the following section, we investigate the possible 
progenitors of such an ultra- long event. 

6.1. Thermal components 

The first hypothesis is that the exceptional duration of 
GRB 111209A is due to a component usually not present 
in normal GRBs. Two super long GRBs (GRB 060218, z 
= 0.033, and GRB 100316D, z = 0.059) were associated 
with a supernova sh ock breakout (jCampana et al.l [20061 : 
IStarling et al.l 1201 ID . These two events are also classi- 
fied as sub-energetic ([Howell fc Cowardll2011[ ). with very 
soft X-ray spectra, and their prompt spectra (including 
the steep decay) shows a thermal emission, which is re- 
sponsible for the long duration. We also do observe a 
thermal component at the start of the XRT observation; 
however this component disappears very soon, and most 
of the prompt phase is free of thermal emission. We can 
t herefore discard this h ypothesis. 

IStarling et al.1 (|2012| ) have found in several GRBs as- 
sociated with a supernova a thermal component. They 
explain it as a hot plasma cocoon associated with the 
jet. However, they observe this component in the steep 
decay phase (where we do not see any thermal emission 
from GRB 111209A using XMM- Newton data) which is 
harder and dimmer than the one we detect at the start of 
the prompt phase. Thus, this model cannot be applied 
to GRB 111209A. 

GRB 101 2 25A, another super-long burst 
ijRacusin et al.l I2010D . shows a strong optical ther- 
mal component. The lack of a redshift measurement 
for this burst made difficult to draw any conclusions on 
the explosion mechanism and two very different models 
were proposed: a minor-b ody tidal disruption by a 
neutron star in our galaxy ([Campana et al.l 1201 It) . and 
the merger of a neutron star wi th a helium star in a 
distant galaxy ([Thone et al.ll201lD . 

Apart from the lack of any strong thermal emission at 
optical (and X-ray) wavelengths, GRB 111209A differs 
from GRB 101225A as it is at least 10 3 times more en- 
ergetic (using the most distant scenario of iThone et al.l 
120111 ). suggesting a different origin for GRB 111209A. 

6.2. Magnetars 

In the recent years, evidence has emerged tha t magne- 
tars c ould account for long GRB properties (e.g. iMetzgerl 
2010). This class of progenitor provides a natural expla- 



nation of extended emissions as long lasting energ y in- 
jection (jDaffOisoiraDIlQiil MMmzEj&jm^MLM ■ In- 
deed, if the typical duration of a long GRB powered by 
a magnetar is of the order of 40-100 seconds, atypical 
values of the magnetic field and spin period could allow 
a duration of ~ 25, 000 seconds. However, such atypical 
values cannot fit both the observed valu es of £, m and 
E v . F or instance, using the formalism of iMetzger et al.1 
(2011), reproducing the observed values of E iso and du- 
ration would lead to E p ~ 180 keV, that is a factor of 3 
lower than the rest frame peak energy of GRB 111209A 
(520 keV). Trying to conciliate these values results in 
non physical solutions (e.g. a neutron star spinning at 
a value exceeding the bre akup value, Pu m = 0.96 ms, 
lLattimer fc Prakash|[200l . 

Another possibility would be to consider the long dura- 
tion of the prompt phase as a plateau phase (i.e. assume 
the true prompt phase has been missed by all instru- 
ments). However, the magnetar model indicates that 
once the plateau phase is finished the decay should be 
mono tonic with a decay index of 1 — 2 (DaH'Oss o et al.l 
|20H . This is in contradiction with our findings from 
the light curve (decay of a = 4.9 ± 0.2 with the presence 
o f a following: platea u). 

iTroia et "aLl (120071) has pr oposed a solution (later ex- 
tended bv lLyons et al.l l2~010) to this problem, with a two- 
stage progenitor: the collapsar would produce a magne- 
tar that would collapse into a black hole a few hundreds 
of seconds later. Under this scenario one could expect 
a double plateau, i.e. a non-monotonic decay after the 
first plateau phase. The first plateau, due to the energy 
injection by the magnetar, would then be the observed 
prompt signal. However, again, one needs atypical val- 
ues of the magnetic field and spin period to explain the 
observed properties of this event, and again this would 
produce inconsistent values of Ei SO and E p , not observed 
in this event. Significantly, the "first plateau" has also 
been observed in the gamma-ray band (it is the prompt 
signal), which is itself unique, making GRB 111209A a 
very peculiar event. 

We thus conclude that a magnetar progenitor for GRB 
11 1209 A is not supported by the data. 

6.3. Tidal Disruption Event 

Swift J1644+57 was at first claimed to be a GRB. It 
is thus possible that GRB 111209A is not a GRB but 
a tidal disruption event. This solution is however not 
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very convincing. Firstly, no host galaxy has been de- 
tected, while one may expect a massive black hole to lie 
in a galaxy bright enough to be detected at small red- 
shift. Secondly, th e light curve should decay as £~ 5 / 3 
(jLodato et al.1 12009). This is however not the case (see 
Sec. H}. In fact, Swift J1644+57 did not show the char- 
acteris tic GRB ligh t curve behavior (so called steep-flat- 
steep, iZhansI [2007; clearly visible in the case of GRB 
111209A . For these reasons we can discard this hypoth- 
esis. 

6.4. Giant stars 

The main difficulty in explaining the nature of 
the progenitor o f GR B 111209A is its duration. 
iWooslev fc Hegerl (|2012D investigated scenarios to pro- 
vide a long duration event with an energy budget typical 
of a long GRB. They proposed four scenarios: 

• i) single super-giant stars with less than 10% solar 
metallicity; 

• ii) super-giant systems in tidally locked binaries; 

• iii) pair-instability in very high mass stars resulting 
in collapse to a black hole; 

• iv) helium stars in tidally locked binaries. 

The scenarios ii) and iv) lead, assuming a typical val- 
ues of 1% Mc 2 for the energy conversion efficiency, to 
luminosities orders of magnitude too faint in comparison 
to the that observed in the case of GRB 111209A . Only a 
100% energy conversion efficiency would marginally rec- 
onciliate these numbers. As for the scenario iii), at a dis- 
tance of z = 0.677, one would expect a bright optical su- 
pernova easily detectable with large aperture telescopes. 
Despite follow-up of GRB 111209A by the Hubble Space 
Telescope and the Very L arge Telesco pe, no evidence of 
a supernova was claimed (|Levanll2012D . The above argu- 
ments leave us with scenario i ), a single supergiant star 
with low metallicity. We note iQuataert &: Kasenl (2012) 
used the same kind of progenitor to discuss the link be- 
tween Swift J1644+57 and a GRB. 

The hypothesized progenitors of long duration GRBs 
are Wolf-Rayet stars (stars with the outer layers ex- 
pelled during stellar evolution). When these layers are 
still present, as in low metallicity super-giant stars with 
weak stellar winds, the stellar envelope may fall-back and 
accretion can fuel the central engine for a much longer 
time. In this scenario, blue super-giant stars can produce 
GRBs with prompt emi ssion lasting about 10 4 seconds 
(|Wooslev fc Heger|[20rl . 

Theoretically, a larger star, such as a red super-giant 
star could also be considered, assuming that the outer 
layers are expelled in order to maintain the free-fall 
time scale to the order of 10 4 seconds. For instance, 
a very weak explosion could eject the convective enve- 
lope of a red supergi ant without producing an observ- 
able supernova, (e.g., iDexteEESase^ HqT^) • However, 
in such a case, this layer of matter is still surround- 
ing the progenitor of the GRB, making the surrounding 
medium density profile much m ore complex. As shown 
by e.g. iChevalier et al. I ()2004D . these complex geome- 
tries should imprint their mark on the afterglow light 
curve, and in fact do (see for instance GRB 050904, 



IGendre e t al. 2007). Here we do not see this complex 
afterglow behavior, supporting the suggestion of a sim- 
ple blue super-giant star. 

If one considers the ratio of normal GRBs to ultra- 
long GRBs, this is shown in the Swift archives to be of 
order 200:1. In comparison, the fraction of blue super- 
iants compare d to Wolf-Rayet stars is typically 10:1 
Eldridee et aTll2008l) - therefore ultra-long GRBs should 
be more frequent than normal long GRBs. At first con- 
sideration, these numbers are difficult to reconcile. How- 
ever, the greatest uncertainty of super-giant collapsars 
is whether matter with high angular momentum will be 
ejected by mass loss or w ill fall into the center of the star 
(|Wooslev fc Heger1 f2012). It is plausible that for most of 
blue supergiant stars there is no fall back of the external 
layers, and thus no long emission. Hence, the observed 
and the expected rates may be reconciled by assuming 
the external layers of most blue super-giants do not fall 
into the stellar core. A rare combination of exception- 
ally large rotation and small mass loss at the end of stel- 
lar evolution could lead to a GRB 111209A like event. 
These properties can be both achieved in case of a low 
metallicity progenitor. This would also help reconcile the 
numbers, as most of the local universe blue super-giant 
stars are not observed with low metallicity. One should 
then consider if it is possible to have these low metallicity 
stars in the local Universe. 

6.5. Low metal protogalaxies in the local Universe 

The most putative hypothesis on the above consider- 
ation is the presence in the local Universe of a galaxy 
with a very low metallicity. In fact, blue supergiants have 
short lifetimes, and the only way of having a low metallic- 
ity blue supergiant is to form it from low-metallicity gas. 
We must thus discuss the possibility of having very low- 
metallicity galaxies in the local universe. If such galax- 
ies should have been common in the early Universe, the 
hierarchical galaxy formation model implies that these 
have merged to form the large r galaxies observed to- 
day (see e.g. iKlypin et alJ H999). However, as pointed 
out by several auth ors, such a mechanism is not effi- 
cient at 100 % (e.g. IKlypin et all 119991: iDiemand et"all 
l2007f l. and galaxies with a very low metallicity should 
still exist i n the local Universe . In fa ct, several of those 
are known (Morale s-Luis et all 1201 1|) . and one object (I 
Zw 1 8) has a measured oxygen abundance of 0.05 solar 
one (|Izotov et al.l H997). These are local objects (with 
z < 0.1) accounting for about .1 % of all the local 
galaxies (Mor ales-Luis et all 120111 ). thus reinforcing the 
idea that GRB 111209A is a very rare event. 

Furthermore, as noted by IMorales-Luis et all (120 11|). 
the l uminosity-metallicity relationship (|Leaueux et al.l 
1979) implies that low-metallicity galaxies cannot be ob- 
served in the more distant Universe. The fact that the 
host galaxy of GRB 111209A has not been resolved by 
the Hubble Space Telescope (jLevanl |20H would then 
support this hypothesis: only the GRB afterglow was vis- 
ible, and GRB 111209A traces the location of a putative 
(very) metal poor galaxy at large distance (z=0.677). At 
this distance, this galaxy would not have been detected 
without the GRB which occurred in it. 

7. CONCLUSION 
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We have presented multi-wavelength observations of 
GRB 111209A . From the prompt emission, we have 
shown that this event is the longest GRB ever detected. 
We have also demonstrated that such a burst would have 
been detected up to z = 1.4 as an ultra- long event. The 
lack of any ultra-long GRBs suggests that such events 
are very rare. 

The spectrum shows a very soft thermal emission dur- 
ing the first thousands of seconds of the X-ray observa- 
tions, not observed at later times. In addition, we have 
detected the presence of a puzzling second component 
seen in X-rays during the fast decay and the following 
plateau phase. This spectral signature is very uncom- 
mon and suggests an unusual nature for the progenitor 
of this burst. 

Explaining these spectro-temporal properties is diffi- 
cult in the standard collapsar model, and we propose 
that this event is representative of a new class of rare 
stellar explosion, that of a blue super giant star with a 
very low metallicity. We stress that the presence of such 
a low-metallicity giant star at this distance of our Galaxy, 
although rare, are not impossible. Thus the progenitor of 
this event could represent the closest link e yer discovered 
with rare population III ste llar explosions (jSuwa & Ioka 
1201 It iNakauchi etaLl 12012ft . One should also note that 
this event would not have been detectable at large red- 
shift with the current instrumentation available. Only 
new ins truments with very large detection areas, such as 
LOFT (jFeroci et al.ll2012t i would be capable of such ob- 
servations. However, at redshift 12, GRB 111209A would 
have a duration of 10 5 seconds (i.e. more than one day) 
making such events very difficult to detect. 

The data of GRB 111209A are very rich, and this pa- 
per has only investigated the most prominent part of it. 
In a second paper, we will present a detailed study of 
the afterglow and the broad band spectrum, thus deduc- 
ing the properties of the emission mechanism and of the 
surrounding medium. 
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